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by 


Donald  R.  Smith— ^ 


ABSTRACT 

The  purpose  of  this  report  is  to  arrive  at  a  better  expression  for 
the  relationships  between  coexisting  phases  of  an  analytical  fluid. 

This  expression  is 


,  1  .  / 3D  _C\  ,  / 3D  _C\  2 

a"  ”al  ^  +  \5B  “  3A/1  +  \5B  '  3  A./  al 


(9AD-5BC)  +  27-  ■  (AD-BC)  - 


54A3B 


50AB' 


9  a: 


where 


+ 


DF 

5AB 


25AB 


2  ( 18AD-5BC) 


ap  r  5 


a =  -  1  ; 

1  Pc 


-^3  -  i  ; 


3  P 


1/ 


Mathematician,  Helium  Research  Center,  Bureau  of  Mines,  Amarillo 
Tex. 


' 

V  Htt/Y.  '■ 
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A  =  Pat  ; 


B  =  t-  p  ; 
2  aaa 


C  =  -  p  +  p  ; 
at  aat 


D  =  -  ( -3p  +  p  ,  , 

6  \  aaa  aaaa/ 


E  =  2  Patt  ; 


2p  -  2p  +  p 


at 


aat  aaat 


G  =  ttt  (  12p  -  4p  +  p  ,  , 

24  \  aaa  aaaa  aaaaa/ 


a  =  -  1 

P 


P 

P 


-  i  ; 


t  =  -  1 


and  p  with  subscripts  represents  the  partial  derivatives  of  p  at 
a  =  0 ,  t  =  0 . 


INTRODUCTION 


2/ 


As  a  result  of  a  journal  article  by  Barieau  (1)—  ,  I  became 


2/ 

—  Underlined  numbers  in  parentheses  refer  to  items  in  the  list  of 


references  at  the  end  of  this  report 


interested  in  the  treatment  of  the  expansion  of  the  chemical  potential 
of  an  analytical  fluid  to  the  third  and  fourth  order. 


' 
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An  analytical  fluid  is  one  in  which  the  chemical  potential  can  be 
expanded  in  a  Taylor  series  about  the  critical  point. 

Landau  and  Lifshitz  (6^)  expanded  to  the  first  order, 


in  a  power  series  in  the  variables  of  v  and  t,  where  v  =  V-V^  and 

t  =  T-T  . 
c 

Edwards  and  Woodbury  (3)  tried  to  extend  Landau  and  Lifshitz'  (6) 

2 

expansion  to  the  second  order.  They  included  a  t  term,  which  should 

3 

have  been  omitted,  and  did  not  include  a  v  term.  Giterman  (4)  pointed 
out  this  error,  and  he  correctly  expanded  the  pressure  in  a  series  in 
v  and  t  to  the  second  order. 

Tisza  and  Chase  (8)  expanded  the  chemical  potential  to  the  first 

order  in  the  variables  of  t  and  p  ,  where  t  =  T-T^  and  p  =  p-p  . 

Mistura  and  Sette  (7)  attempted  to  extend  Tisza  and  Chase's  Q3) 

2 

expansion  to  the  second  order.  They  included  a  t  term,  which  clearly 

/  3 

should  have  been  omitted,  and  did  not  include  a  (p  )  term. 

Griffiths  (5)  observed  that  if  Mistura  and  Sette 's  (7_)  expansion 
exists,  this  precludes  the  possibility  of  any  infinite  singularity  in 
the  specific  heat  at  constant  volume. 

Barieau  (JL)  correctly  treated  the  problem  to  the  second  order  by 
expanding  the  chemical  potential  in  a  power  series  in  the  variables 
of  reduced  density  and  reduced  temperature.  He  pointed  out  that  to 


i 


' 


. 


- 


.  • 


•  •  -  ;■  ' 
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arrive  at  the  correct  expression  for  the  slope  of 


1 


2 


-t 


vs 


t, 


it  is  necessary  to  include  additional  terms  in  the  expansion. 


THIRD-ORDER  EXPANSION 


Let 


(1) 


and 


t  =  - 1  ,  (2) 

c 

P  =  | - 1  >  (3) 

c 


where  p  is  the  density,  T  is  the  temperature,  P  is  the  pressure, 
is  the  critical  density,  is  the  critical  temperature,  and  P^  is 
the  critical  pressure. 

The  third  order  expansion  of  the  chemical  potential  in  the 
variables  of  density  and  temperature  leads  to 
2 

=  =  At  +  Ba2  +  Cat  +  Da3  +  Et2  +  Fa2t  +  Ga4  ,  (4) 

P  \dp/T  1+a \da/ 


where 


P- 


is  the  chemical 


..  n  3/ 
potential  .— 


3/ 


Barieau 


(_1)  had  previously  arrived  at 


this  relationship. 


' 
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The  third  order  expansion  contains  terms  of  the  order  of  a  . 

2  2 
With  t  being  of  the  order  of  a  ,  it  is  not  necessary  to  include  at 

3  5  7 

or  at  terms  as  these  would  be  of  the  order  of  a  and  a  ,  respectively. 

The  constants  A,,  B,  C,  D,  E,  F,  and  G,  in  terms  of  the  equation  of 

state,  are  given  by 


A  =  Pat  5 


B  =  77  P 


or  » 

2  aaa 


C  =  ~Pat  +  Paat  5 


d  =  rl-3p  +  p  /  . 

6\  aaa  aaaa/ 


E 


1 

2  Patt  ’ 


F  =  f(2pat  -  2paat  +  paaat)  ; 


G  =  rrri  12p  -  4p  +  p 

24\  aaa  aaaa  aaaaa 


(5) 

(6) 

(7) 

(8) 

(9) 

(10) 
(ID 


where 


at 


aaa 


(auu)  ; 

\dt3a/ a=o }  t=o 


/.3  \ 


\da3/ 


a=0 ,  t=0 


aat 


I  3  ^ 


aaaa 


\at3a  / 


/  4' 

3a4 

\ 


a=0,  t=0 


a=0,  t=0 


(12) 

(13) 


(14) 


(15) 


, 

' 

' 
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1  A3 

3  P.. 


att 


\St  °al  a-0 ,  t-0 


(16) 


aaat 


and 


aaaaa 


(17) 


a=0 ,  t=0 


(18) 


a=0 ,  t=0 


Equality  of  the  chemical  potential  for  coexisting  phases  requires 


3® „dp  =  °  =  At(a3'al)  +  3  B(val)  +  I  Ct(val 

p,  T 

+ 1  D(a3"ai)  +  Et2(a3'ai)  + 1  Ft(vai 


+  5  G(a3'al 


Equality  of  the  pressure  for  coexisting  phases  requires 

n  n„ 


3/dP 


v3p. 


dp  = 


Pi 


0 


(19) 


(20) 


where  the  subscript  1  refers  to  saturated  gas  and  the  subscript  3  refers 


to  saturated  liquid. 
Equation  20  becomes 

I  At(a3*al)  +  4  B( 


Val)  +  3  Ct(a3~al)  +  5  D(a3'al)  +  2  Et2(a3‘al 


+  |  G(a63-a6i;,  -  0  . 


1/4  4 

+  4  Ft\a3~al 


(21) 


' 

■ 
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Dividing  equations  19  and  21  by  ^a^-a^j,  we  have 


and 


where 


111  2 
At  +  —  BEaa  +  —  CtEa  +  —  DEaaa  +  Et 
3  2  4 

+  ~  FtEaa  +  ^  GEaaaa  =  0  , 


11  1  1  12 

—  A.tSa  +  —  BZaaa  +  ~  CtZaa  +  —  DEaaaa  +  —  Et  Ea 

2  4  3  5  2 


and 

Multiplying  equation  22  by  60  and  equation  23  by  120,  we  have 

60A,t  +  20BEaa  +  30CtEa  +  15DZaaa  +  60Et2 
+  20FtEaa  +  12GEaaaa  =  0 


(22) 


1  1 

+  ~r  FtEaaa  +  —  GEaaaaa  = 
4  6 

0  , 

(23) 

Ea  =  a^ 

+  a^  ; 

(24) 

2 

Eaa  =  a^  + 

2 

a^a^  +  a^  ; 

(25) 

„  3  2 

Eaaa  =  a^  +  a^a^ 

2  ,  3  . 

+  a^  +  at  ; 

(26) 

4  3 

Zaaaa  =  a^  +  a^a^  +  < 

2  2  3  4 

a^a^  +  a^a^  +  a^  ; 

(27) 

^  5  4  ,32 

Eaaaaa  =  a^  +  a^a^  +  a^a^ 

,  2  3  ,  4,5 

a3al  a3al  al  ' 

(28) 

(29) 


and 

2 

60A.tZa  +  30BEaaa  +  40CtEaa  +  24DEaaaa  +  60Et  Ea 
+  30FtEaaa  +  20GEaaaaa  =  0  .  (30) 

2 

Solving  equation  29  for  60Et  ,  we  find 

2 

60Et  =  -  60A.t  -  20BEaa  -  30CtEa  -  15DEaaa 

-  20FtEaa  -  12GEaaaa  .  (31) 


-  - 


10 


Multiplying  equation  31  by  Ea,  we  obtain 

2  2 
60Et  £a  =  -  60 At 2a  -  20BZa£aa  -  30Ct(£a)  -  15D£a£aaa 


-  20FtZaZaa  -  12GEa2aaaa  . 

Substituting  equation  32  into  equation  30  yields 

60A.t£a  +  30B2aaa  +  40Ct£aa  +  24D£aaaa  -  60A.tZa  -  20B2a£aa 
2 

-30Ct(2a)  -  15D2a£aaa  -  20Ft2aEaa  -  12GEaEaaaa 

+  30Ft£aaa  +  20GEaaaaa  =  0  . 

Simplifying , 

t[ 10C (4Eaa  -  3EaEa)  +  10F(3Eaaa  -  2EaEaa)]  +  10B(3Eaaa  -  2EaZaa) 
+  3D(8Eaaaa  -  5EaZaaa)  +  4G(5Eaaaaa  -  3EaEaaaa)  =  0 

From  equations  24,  25,  26,  27,  and  28,  we  find 

2  2 
4Eaa  =  4a^  +  4a^a^  4-  4a^  ; 


2  2 
3EaEa  =  3a3  +  6a  aL  +  3a1  ; 

3  2  2  3 

3Eaaa  =  3a^  +  3a^a^  +  3a.^a^  +  3a  ^  ; 


3  2  2  3 

2Ea£aa  =  2a^  +  4a^a^  +  4a3a^  +  2 a^  ; 


4  3  2  2  34 

8£aaaa  =  8a^  +  8a^a^  +  8a^a^  4-  Sa^a^  4-  8a^  ; 

4  3  2  2  3  4 

bEaEaaa  =  5a^  4-  lOa^a^  +  lOa^a^  4-  lOa^a^  +  5a^  ; 


_  5  r  4  r  3  2  _  2  3  4  n  5 

5Eaaaaa  =  5a^  4  aa^a^  +  5a^a^  +  ba^a^  4  ba^a^  +  5a^  ; 


(32) 


(33) 


(34) 

(35) 

(36) 

(37) 

(38) 

(39) 

(40) 

(41) 


' 


' 
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,  ,  5  ^  .  4  ,32  23  ,  4  5 

and  3£a£aaaa  =  3a^  +  6a^a^  +  6a^a^  +  6a^ai  +  6a^a3  +  3a^ 

From  equations  35  and  36,  we  obtain 


4Eaa  -  3£a£a  =  (a^-a^ 


From  equations  37  and  38,  we  find 


3£aaa  -  2Ea£aa  =  (a^-a^  )  £a 


=  (a3-al)‘ 


From  equations  39  and  40,  we  have 


8£aaaa  -  5£a£aaa  =  (  a^-a^  ^3a^  +  ^a3ai  ^al 


From  equations  41  and  42,  we  obtain 

2 

/  v  /  3  2  2  3' 

5Eaaaaa  =  3£a£aaaa  =  (^a^-a^  )  (2a^  +  3a^a^  +  3a^a^  +  2a^y  . 


Substituting  equations  43,  44,  45,  and  46  into  equation  34,  we  have 


10c(a3-ah  +  10f( 


a3-a  ]  Ea 


+  10B ( a3"a1  )  £a 


+  3D^a3-a1^  +  43^  +  3a1 


+  4G(a3-ai)  (2a3  +  +  2a*)  =  0  . 


Solving  equation  47  for  t  yields 

lOBEa  +  3D^3a3  +  4a^a^  +  3a^  +  4G(^2a3 


t  =  - 


3a3al  +  3a3al  +  2a? 


10 (C  +  FEa) 


Squaring  equation  48,  we  have 


t  = 


lOBEa  +  3D^3a3  +  4a^a3  +  3a^  +  4G^2a3  +  3a^a3  +  3a^a3  +  2a^ 


100  (C  +  FEa)' 


(42) 


(43) 


(44) 


(45) 


(46) 


(47) 


.  (48) 

2 


•  (49) 


.  7 
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Equation  29  may  be  written  as 

lOt (6A.  +  3C£a  +  2F£aa)  +  60Et2  +  20B£aa  +  15D£aaa  +  12G£aaaa  =  0  .  (50) 

2 

Substituting  for  t  and  t  in  equation  50,  we  find  for  the  equation 
relating  a^  and  a^, 


/  2  2\  /  3  2  2  3\ 

10B£a  4-  3D( 3a^+4a^a^+3a^  )  +  4g( 2a^+3a^a^+3a^a^+2a^  ) 


C  +  F£a 


6A+3C£a+2F£aaJ 


+  (60E) 


lOBEa  +  3D^3a^+4a^a^+3a^  'j  +  4G^2a^+3a2a^+3a^a^+2a^ 


2  ,  _  3 


100  (C  +  FEa)' 


+  20B£aa  +  15D£aaa  +  12GZaaaa  =  0 


(51) 


Multiplying  equation  51  by  5 (C  +  FEa)  yields 


-  5(C  +  F£a)(6A  +  3C£a  +  2F£aa) 


2  2\ 
lOBSa  +  3D(3a3  +  43^  +  3a]~J 


/  3  2  2  3 

+  4G(2a3  +  3a3a^  +  3a3a^  +  2a^ 


+  3E 


lOBEa  +  3D^3a2  +  4a3a^  +  3a2 


+  4G(^2a3  +  3a3a^  +  3a3a^  +  2a^ 


3  M2 


+  5(C  +  FEa)' 


20B£aa 


+  15DEaaa  +  12G£aaaa 


=  0 


(52) 


Now,  developing  a3  as  a  power  series  in  a^,  expanded  to  the  fourth 


order,  we  have- 


47 


4/ 


The  details  of  this  development  are  given  in  Appendix  A.. 


' 
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a3  =  -  “li1  + 


3D  _C\  /3D 

5B  “  3A7al  \5B  ‘  3A J  al 


+ 


54A3B 


(9AD-BC) 


,  27D  BCE  ,  DF  G 

+  - - o'(AD-BC)  -  — ^  +  ttt  "  - - ^ 

50ABJ  9A3  5AB  25AB2 


(18AD-5BC) 


>  ■  (53) 


Now,  developing 


as 


5/ 


The 


details  of 


this  development  are  given  in  Appendix  B. 


B  2  (9AD-5BC)  3 

3A.  ai  .2  al 

45  A. 


1 

1350A3B 


150B3E  -  (9AD-5BC) (9AD+5BC)  4 

+  3QAB(9AG-5BF)  al 


(54) 


Similarly,  we  find 


_B  2  (9AD-5BC)  3 

3A  a3  "  .-.2  a3 

45  A. 


1 

1350A3B 


150B3E  -  (9AD-5BC) (9AD+5BC) 
+  30AB(9AG-5BF) 


(55) 


Therefore , 


t 


B  2  (9AD-5BC)  3 

5TT  a  •  ~  o  a . 

3A  1  45A2  1 


1 

1350A.3B 


150B3E  -  (9AD-5BC ) (9AD+5BC ) 
+  30AB(9AG-5BF) 


4 


a . 

l 


(56) 


where  i 


1  or  3 . 


, 

■ 


' 
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Now,  developing  and  a^  as  a  power  series 


1 -  ,  •  6/ 

y  =  ^/-t  ,  we  obtain— 


in  y,  where 


The  details  of  this  development  are  given  in  Appendix  C. 


and 


Now 


al  = 


B 


1  + 


(9AD-5BC)  /-  3At 


30AB 


B 


2  2 

600 A.  B 


(9AD-5BC) (2 1AD-5BC) 
- 2  0 AB ( 9 AG - 5  BF ) 


ioob3e 


3A.t\ 
B  / 


= 


L  3  At 


B 


1  - 


(9AD-5BC) 

30AB 


B 


+ 


2  2 

600A  B 


(9AD-5BC) (2 1AD-5BC ) 
-  20AB(9A.G-5BF) 


ioob3e 


a3+aX 


/-  3  At 

1"  2  ( 9AD-5BC  ) 

F3 At" 

/  B 

30AB 

V  B 

(57) 


(58) 


(59) 


or 


a3+al 


\ 


9AD-5BC1 
5B2 


Thus,  the  reduced  rectilinear  diameter  is  given  by 

1, 


2Va3+a  l.1  +  1  = 


/ 

9AD-5BC 


10B 


t  +  1 


(60) 


(61) 


.  -J.. 


.. 


' 


■ 


■ 
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Then  the  slope  of  the  reduced  rectilinear  diameter  at  the  critical 


point  is  given  by 


L 


d|a3+al 

dt 


9AD-5BC 

10B2 


(62) 


Substituting  equations  5,  6,  7,  and  8  into  equation  62,  we  have 


d( a^+a^ 
”dt 


9p 


at 


t(-3p  +p 
6 \  aaa  aaaa 


-  5(-  p 
\2  raaa 


-pat+P 


aat 


10(l  Paaa)2 


(63) 


or 


d|a3+al 

dt 


5lp  )2 
\  aaa / 


3p  p  -  4p*.  p  -  5p  d 
at  aaaa  at  aaa  aat  aaa 


7/ 


(64) 


7  / 

—  Barieau  (J.)  had  previously  arrived  at  this  relationship. 


From  equations  57  and  58,  we  find 


a 


3 


1 

2  2 

3 00 A.  B 


(9AD-5BC) (2 1AD-5BC ) 
-  20AB(9AG-5BF) 


or 


a 


3 


-a 


1 


Jo  1 

(9AD-5BC ) (2 1AD-5BC )  -  100B  E 

4 

4 

l  3 

100AB 

-  20AB(9AG-5BF ) 

(66) 


Squaring  equation  66,  we 


r 


\ 


obtain 


1 

25AB3 


(9AD-5BC )  (2 1.AD-5BC ) 
-  20AB(9A.G-5BF) 


ioob3e 


(67) 


’ 


..  :  f.  ,  <1  ’  '  J  "■ 
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The  t  term  is  not  included  because  the  coefficient  of  t  obtained  from 
the  fourth-order  expansion  is  different  from  the  coefficient  in  the 
third-order  expansion. 


Dividing  equation  67  by  (-t),  we  have 


Vai 


-t 


3A. 

B 


4  - 


25AB' 


(9AD-5BC) (21AD-5BC) 
-  20AB(9A.G-5BF) 


ioob3e 


(68) 


Thus,  the  slope  of 


a3‘al 


at  the  critical  point  is  given  by 


d_ 

d  1 1 


a3~al 


3  (9AD-5BC) (2 1AD-5BC)  -  100B  E 


-t 


25B 


4  -  20AB(9AG-5BF) 


(69) 


The  second  order  treatment  shows  that  this  slope  is  zero  for  all 
fluids  .  Mistura  and  Sette  (7)  show  that  this  slope  is  constant. 

I  have  carried  out  the  treatment  of  the  expansion  of  the  chemical 
potential  to  the  fourth  order.  That  is, 


2 

fystA  =  _L/3r\  = 

p  \aPy  i+aW/t 


2  3  2  2  1 

At+Ba  +Cat+Da  +Et  +Fa  t 

4  2  3  5 

+Ga  +Hat  -t-Ja  t+Ka 


(70) 


The  coefficient  of  t  in  equation  68  is  the  same  as  the  coefficient 
of  t  obtained  from  the  fourth-order  expansion. 


' 

/ 

, 


‘ 
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For  a  van  der  Waals  fluid,  we  have  (2^) 

P  =  6,  p  =  9,  p=6 

at  aaa  aat 


so  that 


P  =  18,  p  =  0,  p  =  9  , 
aaaa  att  aaat 


and  p  =  45  , 

aaaaa 


9 

A  =  6,  B  =  C  =  0 


D  -  2  s  B  —  0,  F  —  2  > 

and  G  =  . 

o 


Substituting  the  above  values  into  equation  (69),  we  find 


d_ 

dt 


a3~al 


-t 


=8.32  , 


which  agrees  with  the  results  previously  obtained  by  Barieau  (2)  . 


• 

■ 
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APPENDIX  A. 


Let 


X  =  C  +  F£a 


(A- 1) 


Y  =  6A  +  3C£a  +  2F£aa  ; 


(A-2) 


=  10B£a  +  3D^3a^  +  4a^a ^  +  3a ^  +  4G^2a^  +  3a^a^  +  3a^a^  +  2a^  ;  (A-3) 


W  =  20B£aa  +  15D£aaa  +  12G£aaaa  ; 


(A-4) 


or, 


X  —  C  +  F( a^+a^  /  > 


(A-5) 


Y  =  6A  +  3c(a3+aJ  +  2f(^3  +  a3ai  +  a2)  ; 


(A-6) 


Z  -  10B^a3+a^  +  3D^3a3+4a3a^+3a^  +  4G^2a3+3a3a^+3a3a^+2a^^  >  (A-7) 


and 


W  =  20B^a3+a3a1+a^  +  lSD^+a^a^+a^+a^  +  ^G^+a^aj+a^a^+a^  2+a^.  (A-8) 


Let 


dX  _  / 

da  ~  X  ’ 


d  X  =  x„ 


da 


1 


and 


dY 

da. 


=  Y 


d2Y 

da? 


=  Y" 


,  etc 


Now 


(X)a3=0,  ax=0  G  ; 


(A- 9) 


20 


(X  ^=0,  a^O  F\3+1  '  ’ 


(A- 10) 


(X")  _ 


=  Fa! 


a3=0,  a^O  “*3  ’ 


(A - 11) 


(x'")  _ 


=  Fa 


a3=0,  a1=0  ‘“3  ’ 


(A- 12) 


(X  ) 


=  Fa 


//// 


a3=0,  3^=0  ‘"3  ’ 


(A-13) 


(Y)a  =0,  a  =0  =  6 A  ; 


(A- 14) 


(Y\=0,  ai-0“  3Cla3+  l>  ’ 


(A-15) 


(Y’')a3=0,  a^O  3Ca3  +  2F(2a3a3  +  2a3  +  2)  : 


(A- 16) 


(Y")  -  =  3Ca"' +  2F(6a'a"  +  3a'  , 

a3=0,  a^=0  3  V  3  3  3/ 


(A-17) 


/// 


(Y""  )  n  _  =  3Ca""  +  2F  8a  a!  +  6a"a"  +  4a'  , 

a3=0,  a^=0  3  \  3  3  33  3/ 


(A- 18) 


(Z)a3=0,  a^O  0 


(A -19) 


(Z\=0,  ai-0  ■  10B(a3+  0  • 


(A.-20) 


(Z//)a3=0,  ar0=  lOBaJ  +  3D(6a^  +  8a^  +  6  #  , 


(A-2 1) 


(Z'")  _n  _n  =  lOBa!7  +  3D  18a'a"  +  12a" 
a  -0,  3^=0  3  \  3  3  3 


+  4G 


3  2 

12(33)  +  18(as)  +  18a3  +  12 


(A-22 ) 


...  i 


■ 

. 
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/// 


(Z"")  ,  _n  =  10Ba/'//+  3D  24a;a"/  +  ISa^a"  +  16a. 

a^-U ,  a^— U  3  \  3  3  33  3 


+  4G(J2a3a^  +  72a3a"  +  36a3  ,  , 


(W)a3=0,  a^O  =  °  ; 


(W  )a3=0,  a3=0  °  ’ 


(W//)a3-0,  ai=0  =  20Bl2a3a3  +  2a3  +  2^  ’ 


(w  >3^=0,  a^O  2°B(6a3a3  +  3a3  )  +  15D 


6'  [lf+  6(a3)2+  6a3  +  6  * 


and 


(W"")  _n  _n  =  20B 
a3=0,  aL-0 


L8a3a3  +  S  +  Ha3. 


+  15D 


36(a3)2  +  24a^  +  12a;' 


+  12G 


24(a3)4  +  24(a3)3  +  24(as)2 


+  24a3  +  24 


Substituting  equations  A- 1 ,  A-2,  A-3,  and  A-4  into  equation  52,  we 
have 


2  2 

-5XYZ  +  3EZ  +  5X  W  =  0  . 


(A-2 3) 

(A-24) 

(A-25) 

(A-26) 

(A-2 7) 


(A-28 ) 


(A.-29) 


, 


.  .  ■ 
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Thus , 


-SCXYZ'+XZY'+YZX ')  +  6EZZ '  +  5  (XW  7+2XWX7 )  =  0  ;  (A- 30) 


-5 (XYZ77+XZY/7+YZX77+2XY  Z ' +2YX' Z ' +2ZX; Y7 )  +  6E  ZZ77+(Z7) 


/  ,„  / 


r  /  .2 


+5 


X2W/;+4XX '  W '  +2XWX77+2W  (X 7  )  2 


=  0  ; 


(A-31) 


-5  (XYZ///+XZY///+YZX///+3XY 7 Z "+3YX 7  Z/7+3XZ  7  Y77+6X ’  Y'  Z  7+3YZ  7X/7+3ZX '  Y77+3ZY  7X77 ) 


'ttH 


'  ,r  // 


+6E(ZZ///+3Z'Z//)  +  5 


and 


x2w777+2xwx777+6xx  7w77+6xw  7x/7+6w7  (x7  )2+6wx7x 


=  0  ;  (A- 32) 


-5  (X YZ  7777  +XZY7774yZX77/4-4XY 7  Z//7+4YX  'Z///+4XZ/Y///+4YZ  7X777+4ZX  7  Y7/7+4ZY 1  x!" 


+ 1 2X ''  Y '  Z  77+ 1 2X 7  Z 7  Y  77+ 1 2  Y 7  Z 7  X  77+6XY 77  Z  77+  6  YX 77  Z  77+6  ZX 77  Y 77 ) 


/  / 


/  n  !  \7  //_!_  r  \r\7  It  n  //  ,  ^  xrr7  //«//,  r  n\r  Hxt" 


+  6E 


ZZ  7777  +4Z '  Z+3(Z) 


+  5 


X2W ///7  +  8XX  ‘  W 7"+2XWX  7777  +8XW  '  X"+8WX '  X "'+  12XX"W' 


i  ,  2 


+12  (X 7  )^W  77+2 4X  7W 7 X  /7+6W  (X 77 ) 


=  0 


(A- 33) 


Substituting  equations  A-9,  A.-14,  A-19,  A-20,  A.-24,  and  A-25  into 
equation  A-30,  we  obtain 


-5C(6A.)(10B)(a3+l  )  =  0 


(A.-34) 


or 


a3  _1  ‘ 


(A- 35) 


Substituting  equation  A.-35  into  equations  A- 10 ,  A- 15,  A-16,  A, -17, 
A- 18,  A-20,  A-21,  A-22,  A-23,  A-26,  A-27,  and  A-28,  we  have 


“'V”.  v» 


0 


(A- 36) 


h . 

/ 

■ 


- 


■ 


. 
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(Y  >a3=0,  a^=0  0  ; 


(Y">a3=0,  ar0  =  3Ca3  +  4F  > 


(yV°>  ai-o  - 3ca;  - 6Fa;  5 


<Y  >■  -0,  a^O  '  3Ca3  +  4F 


o  //  r\  /// 

3/  ~  3. 


(Z  )a3=0,  a3=0  °  ; 


<Z'\=0,  a1=0  ■  10S  +  120  * 


(Z*>a  -0,  ai-0  “  10Ba3  -  18D«5  ; 


(Z"")  _ 


=  lOBaf  +  6D 


a^=0,  a^=0 


9  2-  4a"' 


+  144Ga"  ; 


(W")  _n  _n  =  40B  ; 

a^-0,  a^-0 


(W  )  _n  _n  =  -  60Ba 

a3-°>  ar° 


3 


and  (W"") 


a^=0 ,  a^=0 


=  40B 


//  \2  0  /// 


3  a3/  -  2a3 


+  360Da"  +  288G 


Substituting  equations  A-9,  A-14,  A-19,  A-20,  A-24,  A-36,  A-42, 


(A- 37) 

(A- 38) 

(A-39) 

(A- 40) 

(A-41) 

(A-42) 

(A-43) 

(A -44) 

(A- 45) 

(A -46) 

(A.-47 ) 

A -45 , 


and  A.-31,  we  obtain 


24 


-5C(6A)(lOBa"  +  12d)  +  5C2(40B)  =  0 


(A-48) 


or 


//  =  of—  — 

L3  \5B  3  A. 


(A-49) 


Similarly,  equation  A.-32  becomes 


-5C(6A)(lOBa"'  -  18Da")  +  5C2(-60Ba")  =  0  ; 


(A-50) 


thus 


a-  =  (m  .  c\  „ 

3  V5B  A./  3 


(A-51) 


or 


■  -  «8§  -  £)' 


(A.-52) 


Similarly,  equation  A-33  becomes 


■5C(6A)|lOBa""+  6D 


9(  a3  J  ‘S' 


+  144Ga 


"}  -  30c(3Ca"+4F^10Ba"+12D) 


-30(6A.)(Fa^(^10Ba"+12D^  +  18E  (lOBa"+12D)2+  5C2|40B 


3(  a3/  ~  2< 


+  360Da"  +  288G  j-  +  60C^Fa"^(40B)  =  0  . 


(A-53) 


Equation  A-53  may  be  written  as 

75  ABC  a^  -  2C>(9ACD-5BC2)a"/  +  15 (30ABF+27ACD-30B2E+5BC2)(a")2 

+  60^18ACG+9ADF-5BGF-18BDE-3C2D)a"  +  72 (5CDF-5C2G-9D2e)  =  0  .  (A-54) 


<■ 

!  ; 


.  .•  ■ 1 • : •  - 
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Substituting  equations  A-49  and  A-52  into  equation  A-54  yields 

25 ABCa/'//+  -8C  C9AD:5^C  )-  +  — —  "^C).  (30ABF+27ACD-30B2E+5BC2) 

3  45A  B  45A  B 

“  8"9t^'5BC  ( 18ACG3-9ADF-5BCF-  18BDE-3C2D)  +  24(5CDF-5C2G-9D2E)  =  0.  (A-55) 

jAB 

Simplifying, 


lOOC' 


324CD 


25ABCa"//  +  (9AD-BC )  +  (AD-BC) 


9  A. 


B 


2  2 

200B  C  E 
3A2 


+  120CDF  -  —  (18AD-5BC)  =  0  . 

B 


(A-56) 


Solving  for  a "" ,  we  have 


////  — 


4C 

3 

9  A,  B 


(9 AD-BC )  - 


324D  (AD_fiC)  +  8BCE 


2  5AB~ 


3a; 


24DF  +  -l4^  (18AD_5BC) 


5AB 


25AB 


(A.-57) 


We  now  have 


the  equation  for 


as 


a  function  of 


expanded 


to 


the 


fourth  order, 


a 


3 


3A;ai  + 


/3D 

\5B 


a 


2 

1 


+ 


54A3B 


(9 AD-BC)  + 


27D  ,  _  BCE  ,  DF 

- -  (AD-BC)  "  — 7  +  77T 

50AB3  9  A,3  5AB 


-  — —■  (  18AD-5BC  ) 
25AB 


(A.-58) 


,  '  • 
'  \‘v  > 


- 
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Let 


and 


APPENDIX  B 


=  3a2  +  4a^a^  +  3a2 


3  2  2  3 

e  =  2a^  +  3a^a^  +  3a^a^  +  2a^ 


(  B  —  I ) 


(B-2  ) 


Substituting  equations  24,  B-l,  and  B-2  into  equation  48,  we  have 


Let 


Now 


t  =  - 


lOB^a  4-a^  +  3D6  +  4Ge 


10 


C  +  F( a3+a^y  , 


^3  _  / 

da3  a3  ’ 


d2a, 

3  _  // 

2  a3  ’ 
da^ 


d6  


da. 


=  6 


da2 


a7', 


.  ,  etc 


2 (9AD-5BC) 
15  A  B 


2 (9AD-5BC )2 

2  2 
75 A  B 


(B-3) 


(B-4) 


(B-5 ) 


(B-6) 


27 


//// 


=  -  24 


a^-0 


C  (9AD-BC)  +  ■  27D0  (AD-BC)  -  +  °F 


54A3B 


50AB" 


9  A 


3  5AB 


25AB 


( 18AD-5BC ) 


(6)a3=0,  a^O  °  5 


(6/)a3=°,  ax=0  0  5 


^a^O,  a^O  6 \a3/  +  8a3  +  6  5 


<*' Vo,  ai=0  =  18a3a3  +  12a3  > 


(6  )  _  _ 


a3=0,  ai=0  =  2Sa;  +  18W;  +1S  ’ 


(e)  _n  _n  ~  0 

a3-°,  ar0 


(e  4^0,  a^o  0  ; 


(e")a3=0,  ax-0  =  °  > 


3  2 

(e"')l  =Q  a  „  =  12(ai)  +  18(4:4  +  18al  +  12 


a  =0,  a  =0  “V3)  '  "V~3/  '  ‘”“3  '  -  ’ 


(B-7) 

(B-8) 

(B-9) 

(B- 10) 

(B-ll) 

(B- 12 ) 

(B-13) 

(B- 14) 

(B- 15) 

(B-16) 


- 


/ 
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and 


<e/'")a3=0>  a1=0  ■  72K)  a5  +  72a3a3  +  3S  ’ 


(B-17) 


Substituting  equations  B-4,  B-5 ,  and  B-6  into  equations  B-10,  B-ll, 
B-12,  B-16,  and  B-17,  we  have 


(6'V°’  ai=o  =  4 


(6'")  _ 


a^=0 ,  a^=0 


=  4(9AD-5BC) 

5AB  ; 


(B- 18) 


(B-19) 


(6"")  _ 


a3=0,  a^O 


8(9AD-5BC)‘ 

2  2 
15A.  B 


(B-20) 


(* >a3-0,  ai-0  “  °  J 


(B-21) 


and 


(s  ).  - 


a3=°’  al“° 


24(9AD-5BC) 

5AB 


(B-22) 


From  equation  B-3,  we  have 


t  =  - 


10B 

a3+l^  +  3D6  +  4Ge 

lOB^+a-J  +  3D6  +  4Ge 

F 

Va3+1) 

10 

"C  +  F(a3+al)” 

I 

10 

>  +  F(a3+al)! 

2 

5  (B-2 3) 


10Ba"  +  3D§ "  +  4Ge"  2  lOB^a'+ly  +  3D6 '  +  4Ge'  F^a'+l) 


>£c  +  F^ 


10|  C  +  F(  a3+a1 


—  + 


10 


C  +  F(a3+ai 


+ 


10B(a3+a1j+  3D6  +  4Ge](Fa"  )  2 


10B^a3+ai j  +  3D6  +  4Ge 


F( a3+! 


10 


C  +  F( a3+ai 


2 


10 


C  +  F(a3+ai 


_3 


(B-24) 


•  r 


' 


,  •  . 
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t  =  - 


lOBa"'  +  3D  6  +  4  Ge 


/// 


10 


C  +  F( a^+a^ 


+ 


3  lOBa"  +  3D§ "  +  4Ge 


F  a„+l 


10 


C  +  F( 


+ 


10B^a^+l J  +  3D6  +  4Ge 


Fa: 


10B 


(a3+0 


_2 


+  3D6  +  4Ge 


F(a'+1 


10 


C  +  f( a.+a,  ) 

10 

L  V  3  1/J 

C  +  F( a^+a^  , 


+ 


10B(^a3+a1)  +  3D6  +  4Ge 


Fa: 


10B(^a3+a1  )  +  3D6  +  4Ge 


F(a3+1 


Fa; 


10 


C  +  F( a3+aL 


10 


C  +  F(a3+a 


+ 


10B^a3+a]L  )  +  3D6  +  4Ge 


F'a3+1 


-.3 


10 


C  +  F(a3+a1 


_4 


and 


(B-25) 


t  ""  — 


lOBa^H-  3D6///;+  4Ge""  4(jL0Ba3  +  3D6'"  +  4Ge'")  F(a'+l) 
_  +  —  2~" 


10 


C  +  F(  a3+a^ 


10 


C  +  F( a2+a^ 


+ 


6(^10Ba3  +  3D§ "  +  4Ge"J(^Fa3j  12^10Ba"  +  3D6  "  +  4Ge 


Ffa'+l 


10 


C  +  F(  a3+a;L 


10 


C  +  F( a3+a^ 


24 


+ 


lOB^a'+l )  +  3D6  '+  4Ge 


F(a3+1 


ic 


10|  C  +  F(  a3+a3 


_4 


+ 


10B(^a3+a1)+  3D6  +  4Ge 


Fa 


//// 


10 


r 


C  +  F( a3+aL 


8 


10B(a3+a1  )  +  3D6  +  4Ge 


F(a3+1 


w 


Fa3  1  6 


10B(a3+a1  )+  3D6  +4Ge 


Fa: 


10 


C  +  Fl a3+a^ 


_3 


10 


C  +  F( a3+a^  j 


-,3 


36 


+ 


lOB  (  a0+a1  )  +3D6  +4GeTF(a'+l)  ^(Fa"  )  24 


\  3  1 


JL  V  3 


10b ( a3+a1 


)+3D6+4Ge  Jj^a'+l 


|4 


10 


C  +  F( a3+a3 


4 


10 


c  +  F( a3+a1y 


-5 


4 


lOB  ^a^+l^+3D§  +4Ge 


/// 


Fa"'  )  24 


10B(a3+lj+3D6/+4Ge/ 


F  a3+1 


Fa 


10 


C  +  F( a3+a| 


_2 


10 


C  +  F(a3+a1 


-,3 


.  (B-26) 
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Substituting  equations  B-8  and  B-13  into  equation  B-3,  we  obtain 


(t) 


v°’ 


a^=0 


=  0 


(B-27) 


Substituting  equations  B-4,  B-9,  and  B-14  into  B-23,  we  have 


(t') 


al=0 


0 


(B-28) 


Similarly,  equation  B-24  becomes 


(t") 


a3=0, 


10B 

r-2(9AD-5BC)n 

15AB 

+  3D(4) 

-10C 

(B-29) 


or 


V° 


2B 
3  A  * 


(B-30) 


Similarly,  equation  B-25  becomes 


ar° 


10B 


2 ( 9AD-5BC) 

2  2 
75A  B 


+  3D 


~4(9AD-5BC ) 

5AB 


-10C 


(B-31) 


or 


(t") 


a^=0 


2 (9AD-5BC ) 
15  A2 


(B-32) 


. 


' 


1  .  .  1 


■  . 
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Similarly,  equation  B-26  becomes 


(t"")  _ 


a3=0,  a1=0 


10B(-24) 

-IOC 


2  2 
C  v  27D  ,  x  BCE 

(9AD-BC)  +  — — r  (AD-BC)  -  — ~ 


54A3B 


50AB' 


9A' 


+ 


DF 

5AB 


25AB 


2  ( 18AD-5BC ) 


3D 


8  (9AD-5BC  )‘ 

2  2 
15A  B 


+  4G 


-24(9AD-5BC) 

5AB 


+ 


-IOC 


or 


6F 


+ 


-2 (9AD-5BC ) 

15  AB 


iob(- 


-2 (9AD-5BC ) 

15  AB 


+  3D (4) 


IOC' 


,  (B-33) 


(t"")  _ 


4 


(9AD-5BC  )  (9AD+5BC  ) 


a3=0,  a^O  225A3B  L"150BjE  -  30AB(9A,G-5BF). 


(B-34) 


We  now  have  the  equation  for  t  as  a  function  of  a^, 


-B  2  (9AD-5BC )  3 

t  =  —  a,  -  - ~ — —  a 

3A.  1  /ca2  1 

45  A 


1350A.  B 


150B3E  -  (9AD-5BC ) (9AD+5BC ) 
+  30AB(9AG-5BF ) 


at.  (B-35) 


. 


■/ 


I 


I 
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Let 


APPENDIX  C 


Let 


a  = 


-A 

3  A.  5 


(C-l) 


3  = 


9AD-5BC 
45  A  2 


(C-2  ) 


and 


Y 


1350 A  B  L. 


150B  E  -  (9AD-5BC) (9AD+5BC) 
+  30AB(9A.G-5BF ) 


(C-3) 


Thus , 


2  3  4 

-t  =  cva.  +  Ba.  +  va. 

l  l  '  i 


(C-4) 


-t  =  y 


(C-5) 


or 


y  =  ±a/- t  ,  (t  <  0)— ^  , 


(C  -  6 ) 


— !  t  is  below  the  critical  point  and  is,  therefore,  negative 


where  the  minus  sign  in  front  of  the  square  root  applies  to  the  saturated 
gas  (a^)  and  the  positive  sign  in  front  of  the  square  root  applies  to 
the  saturated  liquid  (a^) . 


Differentiating  equation  C-4,  we  find 


da . 


dt  9  2  3 

2o,a.  +  3Ba.  +  4va. 

l  i  1  i 


(C-7) 


/ 


; 
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Differentiating  equation  C-5,  we  find 


-  ^  _  2y  -  2{aal  +  ^ai  +  yal) 


1/2 


Thus , 


da  3 


2\ 


1/2 


dy 


i\  _  2\?  +  Pai + 


2(y  +  3|3a .  +  4ya . 

l  1  l 


and 


Id2  \ 

d  a . 

l 


,  2 
\dy  / 


^@+2yai^  4^cH-(3ai+ya2^33+8ya_^ 

^2o+3@a^+4ya2^2 


(^2o+33ai+4ya^^ 


ld\] 

8y(af3a.+ya3) 

1/2  16 

^3+2yai^3g+8yai^ 

(o+ga^ya3)17'2 

\dy3  J 

( 2o+3|3a .+4ya3 

3 

^2o+33a  .+4ya2^ 

^2^3/2 
i  '  i. 


^2o+3f3ai+4ya^ 


4 


Now 


'da  .> 

l) 


^'a.=0  > 


U2,\ 

i 


\dy  /a  = 


a .  =0 
i 


£_ 

2  5 
a 


and 


/  j3 

d  a . 

l 


\dy 


+ 


15£ 


a .  -0 

l 


5/2  7/2 

o'  4cy 


(C-8) 


(C-9) 


(C-10) 


(C-ll) 


(C-12) 


(C- 13) 


(C- 14) 


, 


' 
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We  now  have  the  equation  for 


a . 

l 


as 


a  function  of  y, 


or 


+ 


/  7/2, 

4a 


3 

y 


(C-15) 


(C- 16) 


For  i  =  1,  y=-/v/-t,  so  that 


(C-17) 


Substituting  equations  C-l,  C-2,  C-3  into  equation  C-17,  we  obtain 


a 


1 


(9AD-5BC) 

30  A  B 


+ 


1 

2  2 

600A.  B 


(9AD-5BC)(2 1 AD - 5  BC ) 
-20AB(9AG-5BF) 


ioob3e 


-3A.t\ 


B 


; 


(C- 18) 


y 


For  i  =  3 ,  y  =  +  //-t  ,  so  that 


if  -Ft  -  if  (v  -  Hr) 


(C- 19 ) 


Substituting  equations  C-l,  C-2,  and  C-3  into  equation  C-19,  we  have 


a 


3 


(9AD-5BC )  j-JAt 

30AB  V  B 


,  1 

(9AD-5BC) (2 1AD-5BC )  -  100B3E 

f-3At\  , 

1  2  2 
600A.  B 

-  20AB(9AG-5BF) 

\  B  / 

- 

(C-20) 


